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Abstract 

We demonstrate that the decay r] Sn represents a sensitive probe 
for the breaking of chiral symmetry by the quark masses. The transi- 
tion amplitude is proportional to the mass ratio (m^ — m^) / (m^— m^). 
The factor of proportionality is calculated by means of dispersion re- 
lations, using chiral perturbation theory to determine the subtraction 
constants. The theoretical uncertainties in the result are shown to be 
remarkably small, so that ?7-decay may be used to accurately measure 
this ratio of quark masses. 
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1. Chiral perturbation theory. Since Bose statistics does not allow three 
pions to form a configuration where the total angular momentum and the 
total isospin both vanish, the decay 77 — Svr violates isospin symmetry. The 
transition amplitude contains terms proportional to the quark mass difference 
^d~fT^u as well contributions oc e^, generated by the electromagnetic interac- 
tion. The latter are strongly suppressed by chiral symmetry - the transition 
is due almost exclusively to the isospin breaking part |(m„ — mrf)(uM — dd) of 
the QCD Hamiltonian Current algebra implies that the matrix element 
relevant for the transition rj tt+tt^tt'^ is given by 



^ = -^0^(.-3Mj , ^0 = -^-3-y, (1) 



where = 92.4 MeV is the pion decay constant, s = + Pn-Y and rh 
denotes the mean mass of the u- and (i-quarks, rh = |(mu + m^). 

The current algebra formula represents the leading contribution of the 
chiral perturbation series, which describes the decay amplitude as a sequence 
of terms involving increasing powers of momenta and quark masses. The 
corrections of first nonleading order (chiral perturbation theory to one loop) 
are also known M. If the amplitude is written in the form 



1 M|M|.-M2 1 _ml-ml 



the one loop result for the dimensionless factor M{s, t, u) exclusively involves 
measured quantities. It is of the form 

M(s,t,u) =Mo(s) + (s--u)Mi(t) + (s-t)Mi(M)+M2(t)+M2(M)-|M2(s) , (3) 

with t = {pt,- +Pt,oY, u = {pt,o +p^+)^. The functions Mo(s), Mi(s), M2(s) 
represent contributions with isospin I = 0, 1, 2, respectively. They contain 
branch point singularities generated by the final state interaction in the S- 
and P-waves. The imaginary parts due to partial waves with angular momen- 
tum i > 2 only start showing up if the chiral perturbation series is evaluated 
to three or more loops. For this reason, the amplitude retains the structure 
@ also at next-to-next-to-leading order. We stick to this approximation 
throughout the following, i.e. represent the amplitude in terms of the S- and 
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P- waves and neglect the discontinuities from partial waves with £>2.Q 

In current algebra approximation, only the 1=0 component Mo(s) is dif- 
ferent from zero. According to eqs.(l) and (2), the leading term is given 
by T{s) = (3s - 4M2)(M2 - M^)'^. The one loop calculation supple- 
ments this expression with corrections of first nonleading order and also gives 
rise to components with 1=1,2. We denote the one loop approximation by 
Mo(s), Mi(s), M2(s). In the notation of ref.p, the result reads 

Mo{s) = T{s) + Ao(s){l + |T(s)} + A2(s){l - |T(s)} + A,{s) + V{s) 

+ |AGMor(s) + {ApT(s) - Agmo}(M2 - M^)-' , 
Mi(s) = {|Ai(s) - AsL,F-'}iM'^ - Mir' , 
M2(s) = iA2(s){3-T(s)} . 

The functions Ao(s), . . . , A3(s) are generated by the final state interaction. 
They arise from one loop graphs, contain cuts along the real axis and are 
inversely proportional to F^. The quantity V{s) represents a subtraction 
term linear in s, which also originates in these graphs (for explicit expressions, 
see ref. 0). The coefficients Aqmo and Ap are known low energy constants 
related to the masses and decay constants of the pseudoscalar octet. Finally, 
the term L3 is an effective coupling constant of the chiral Lagrangian, whose 
value is known from fr^4-decay 0. 

The one loop formula shows that, up to and including first order correc- 
tions, the decay rate is inversely proportional to Q^, with a known factor of 
proportionality. A measurement of the rate thus amounts to a measurement 
of the quark mass ratio Q. Apart from experimental uncertainties, the ac- 
curacy of the result for Q is determined by the accuracy of the theoretical 
prediction for the factor M{s,t,u). This motivates the present work. 

In fig.l, the real part of the one loop amplitude is plotted along the line 
s = u (dash-dotted). For small values of s, the curve very closely follows 
the current algebra formula (dashed) - there, the correction is small. The 

^The same approximation is commonly used also in tttt scattering, where the Roy 
equations |^ provide a rigorous starting point for the dispersive analysis. The validity of 
the representation analogous to (|^) to two loop order of chiral perturbation theory was 
pointed out in ref.Q and an analysis of the tttt data on this basis is given in ref.Q. The 
explicit two-loop result not only confirms this structure, but shows that reasonable 
estimates for the relevant effective coupling constants lead to very sharp predictions for 
the scattering lengths and effective ranges. 
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Figure 1: Real part of the amplitude along the line s = u. The dashed line 
is the current algebra prediction. The dash-dotted curve corresponds to the 
one loop approximation of chiral perturbation theory. The full line represents 
the result of the dispersive calculation described in the text. 

cusp at s = 4 arises from the final state interaction, which is dominated 
by the isospin zero S-wave. The singularity requires the amplitude to grow 
more strongly with s than the current algebra result. At the maximal energy 
accessible in ?7-decay, ^/s = — M^, the correction amounts to | of the 
leading term. This indicates that (a) the current algebra formula underesti- 
mates the rate by a substantial factor and (b) the first two terms of the chiral 
perturbation series represent a good approximation only at small values of s. 

The reason why in part of the physical region, the corrections are un- 
usually large is understood. Chiral symmetry implies that the pions are 
subject to an interaction with a strength determined by F^. In the chiral 
limit the corresponding 1=0 S-wave phase shift is given by 6q = s/lGvr^F^ 
and thus grows with the square of the center of mass energy. Accordingly, 
the corrections rapidly grow with the energy of the charged pion pair. 

2. Dispersion relations and sum rules. There is a well-known method 
which allows one to calculate the final state interaction effects, also if they 
are large: dispersion relations. Indeed, this method was applied to three- 
body decays long ago and many papers concerning the decay 77 Stt have 
appeared in the meantime - they may be traced with refs.[§-||l2|. Crudely 
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speaking, analyticity and unitarity determine the decay amplitude up to the 
subtraction constants. Chiral perturbation theory is needed only for these. 
The additional information contained in the one loop result merely reflects 
the fact that the chiral perturbation series obeys unitarity, order by order. 
In the following, we determine the absorptive part with unitarity and use 
dispersion relations to express the full amplitude in terms of the absorptive 
part and a set of subtraction constants. The latter are calculated by matching 
the dispersive and one loop representations at low energies. In principle, the 
same steps must be taken whenever unitarity is used to extend the range of 
validity of the effective theory - our analysis outlines the general procedure 
in the context of a nontrivial example. For a recent report on the marriage 
of chiral perturbation theory and dispersion relations, see ref. 

We first set up the dispersion relations obeyed by the one loop amplitude. 
The functions Mo(s), Mi(s), M2{s) are analytic except for a cut along the 
real axis. For s —* oo, Mq{s) and M2(s) grow like s^logs, while Mi(s) is 
of order slogs. Accordingly, three subtractions are needed for Mq{s), M2(s) 
and two for Mi(s), 

— , , , r, r°° ds' ImMofs' 
Mo s = ao + 60S + cos^ + - / —- ^ 

TT J4Af2 S'-^ S' — S — 2 



— , ^ , s^ /-oo ds' ImMi(s') 

TT J4M| S'^ S' - S - le 

2 s^ ("^ ds' ImM2(s') 



M2(s) = 02 + 62S + C2S^ + — f 

IT J4 



vr JiM^ s'^ s' — s — ie 

Inserting this representation in eq.(^), the one loop amplitude takes the form 
M(s, t,^) = P{s,t,u) + D{s,t,u), where the polynomial P{s,t,u) contains 
the contributions from the subtraction constants ao, . . . ,C2 while D{s,t,u) 
collects the dispersion integrals. On kinematic grounds, the polynomial 
only involves four independent terms, which may be written in the form 
P(s, t,u) = a + b s + c s"^ — d {s'^ + 2 1 u): Only four combinations of the eight 
subtraction constants ao, . . . , C2 are of physical significance. 

The dispersive representation of the one loop amplitude thus appears 
to involve four subtraction constants. In fact, however, the constant c is 
determined by the imaginary parts through a sum rule, for the following 
reason. The amplitudes Mo(s) and M2(s) individually grow like s^ logs for 
s — *■ 00, but the combination Mo(s) + | M2(s) only diverges in proportion 
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to slogs. For the above representation to reproduce this behaviour at large 
values of s, the term oc s^ arising from the dispersion integrals must be 
compensated by the one from the subtraction polynomials, i.e. 

1 ds' 

c = Co + |c2 = - / -r{ImMo(s') + I ImM2(s')} . (5) 

The terms of order s^ logs also cancel in the combination s Mi(s) + M2(s), but 
this function does not obey a twice subtracted dispersion relation, because 
the term oc L3 gives rise to a contribution which grows like s^. Hence d is 
related to the effective coupling constant L3, 

^ - ^.+c. ^ - '^^y^g +^ ^{.■ImM.(.-)+ImM.(.-)} . (6) 

This shows that the one loop amplitude is fully determined by its imag- 
inary parts and by the three constants a, h and L3. The dispersive repre- 
sentation explicitly expresses the fact that an analytic function is uniquely 
determined by its singularities: While the functions ImMi(s) describe those 
associated with the branch cut along the real axis, the subtraction constants 
account for the singularities at infinity. 

3. Unitarity. In one loop approximation, the amplitude obeys unitarity 
only modulo contributions of order p^. The corresponding expressions for 
the discontinuities across the branch cut only hold to leading order of the 
chiral expansion. The unitarity condition may be solved more accurately, 
using an approximation which holds up to and including two loops of chiral 
perturbation theory. The approximation accounts for two-body collisions, 
but disregards multiparticle interactions. In the case of elastic scattering, 
the approximation is well-known and is referred to as elastic unitarity. The 
extension of this notion to three-body decays is not trivial, however: Some 
of the early papers contain erroneous prescriptions. A coherent framework 
is obtained with analytic continuation in the mass of the rj. One first con- 
siders the scattering amplitude Tit] — > tttt for a value of the strange quark 
mass for which Mjj < SM^^, such that the decay is kinematically forbidden. 
Elastic unitarity then represents a consistent approximation scheme. In this 
approximation, only two-pion intermediate states are taken into account, 
such that the imaginary part of the scattering amplitude coincides with 
the discontinuity across the two-pion cut and is determined by the phase 
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shifts of vrvr scattering. Denoting the discontinuity of the function Mi{s) by 
discMi(s) = {Mi(s + ie) — Mi(s — ie)}/2i, the elastic unitarity conditions 
readQ 

discMi(s) = sin 5i{s) e-'^'^'^{Mi{s) + Mi{s)} . (7) 

The first term in the curly brackets accounts for repeated collisions in the 
s-channel, while the second arises from two-particle final state interactions 
in the t- and w-channels and involves angular averages of the type 

«:={(!- 4M^Js){{M, + M^f - s}{(M, - M^f - s}f , 

with so = + M^. In this notation, the explicit expression for the discon- 
tinuity due to two-body collisions in the crossed channels reads 

Mo = |(Mo) + f (M2) + 2{s - so)(Mi) + |/€ (zM,) 

Ml = {3{zMo) - 5{zM2) + |(s - So){zMi) + f (;2'Mi)} 

M2 = (Mo) + |(M2)-|(s-so)(Mi)-i/€(^Mi) . 

In contrast to the imaginary parts, the discontinuities are analytic in M^, so 
that the above relations may be continued to the physical value of the mass. 
While the imaginary parts are real by definition, the two particle disconti- 
nuities obtained through analytic continuation are complex. In the unitarity 
relation for the decay amplitude, these contributions arise from the discon- 
nected part of the T-matrix, describing two-particle collisions with the third 
pion as a spectator. The imaginary part of the amplitude necessarily also 
receives contributions from the connected part. In the above approximation, 
these arise from successive collisions among two different pairs of pions. Note 
also that the dispersion integrals extend beyond the physical region of the 
decay, where k may develop an imaginary part. The integration over z must 
then be deformed into the complex plane, in such a way that the path avoids 
the branch cut (for details, see ref. |Tl|). 

^ In the standard notation, where the phase shift of the partial wave with isospin and 
angular momentum quantum numbers 1,£ is denoted by (5](s), the unitarity condition 
for the amplitudes Mo{s),Mi{s) and A/2(s) involves do{s) = Sq{s),Si{s) = Sl{s) and 
S2{s) = (5o(s), respectively. 
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The conditions (|^) specify the approximate form of unitarity used in our 
analysis. We briefly comment on the three phase shifts 6o{s) , 6i{s) , 62(3) 
which occur in these conditions and which represent an important input of 
our calculation. Chiral perturbation theory yields remarkably sharp predic- 
tions near threshold. Together with the experimental information available 
at higher energies, the Roy equations then determine their behaviour up to 
KK threshold, to within small uncertainties |14]. The function 6o{s) slowly 
passes through ^tt at ^/s = 850 — 900 MeV and continues rising beyond this 
point. In the context of our calculation, the behaviour at higher energies does 
not play a crucial role. We use a parametrization where So{s) tends to vr for 
s — 00, so that the discontinuity, which is proportional to sin^ol^), tends to 
zero. A similar representation of the high energy behaviour is also used for 
the P-wave, where the phase shift rapidly passes through ^vr at s = M^. The 
exotic partial wave with 1=2, on the other hand, describes a repulsive inter- 
action and does not exhibit resonance behaviour. We use a parametrization 
where the function (52 (s) asymptotically tends to zero. 

4- Ambiguities in the solution of the dispersion relations. If the phase 
shifts are known, the unitarity conditions determine the discontinuity across 
the branch cut in terms of the amplitude itself. Analyticity then implies that 
the functions Mq{s), Mi(s) and M2(s) obey dispersion relations analogous to 
eq-d), 

M,{s) = P,{s) + - / ds' si^^il^Jg {M,{s') + M,{s')} . (8) 

TT J4M2 S' — S —te 

For simplicity, we assume that the phase shifts rapidly reach their asymptotic 
values, so that the dispersion integrals converge without subtractions, but 
this is not essential - we might just as well write the integrals in subtracted 
form, absorbing the difference in the polynomials -Pi(s). 

It is crucial that the dispersion relations used uniquely determine the 
amplitude. In fact for the above relations this is not the case: The cor- 
responding homogeneous equations - obtained by dropping the polynomi- 
als Pi{s) - possess nontrivial solutions. In its simplest form, the problem 
shows up if the contributions to the discontinuity from the angular averages 
over the crossed channels are dropped. Unitarity then reduces to three in- 
dependent constraints of the form discMj = sin(5i e~*'^^Mi, or, equivalently, 
Mi(s + ie) =e^*'^^Mi(s — ie). This condition is well-known from the dispersive 



7 



analysis of form factors and can be solved explicitly: The Omnes function, 
defined by 

[n Jam^ s' [s' - s - le) J 

obeys Qi{s + ie) — e^"^' — ie), so that the ratio mi{s) — Mi(s)/Qi{s) 
is continuous across the cut. Since Oi(s) does not have any zeros, mi{s) 
is an entire function. In view of the asymptotic behaviour of the phase 
shifts specified above, no(s),^i(-s) tend to zero in inverse proportion to s, 
while ^2{s) approaches a constant. Assuming that Mi{s) does not grow 
faster than a power of s, the same then also holds for mi(s). Being entire, 
mo{s),mi{s) and ni2{s) thus represent polynomials: The general solution of 
the simplified unitarity conditions is of the form Mi{s)=mi{s) Qi{s), where 
mi(s) is a polynomial. The ambiguity mentioned above arises because the 
Omnes factors belonging to the partial waves with 1=0,1 tend to zero for 
s — > 00. The dispersion relation for Mo(s), e.g. implies Mo(s) — > Po{s). If 
Po{s) is a polynomial of degree n, mo{s) is of dcgcc n + 1. Hence the above 
dispersion relation for Mo{s) admits a one parameter family of solutions. The 
same is true of Mi(s), while the solution of the dispersion relation for M2(s) 
is unique. 

The same problem also occurs if the angular averages are not discarded. 

The preceding discussion points the way towards a solution of the problem: 
It suffices to replace the above integral equations with the dispersion rela- 
tions obeyed by the functions mi(s) = Mi{s) / fli{s) . Since the corresponding 
discontinuities are given by 

Mi{s + ie) Mi{s - ie) _ Mi{s + ie)e-'^' -Mi{s - ie)e'^' _ ^. sin^iMj 
fli{s + ie) Qi{s ~ ie) \Qi\ ' 

the dispersion relations take the form 

X ^ / X ff^/ X 1 /'"^ , / sin 5i(s') Mi(s') ] 

In the simplified situation considered above, where the angular averages 
Mi{s) are discarded, this form of the dispersion relations indeed unambigu- 
ously fixes the solution in terms of the polynomials -Pi(s). Our numerical 
results indicate that the same is true also for the full set of coupled integral 
equations, but we do not have an analytic proof of this statement. 
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In the present context, the asymptotic behaviour of the phase shifts is 
an academic issue. The fact that the solution of eq.(|^) is unique if all of 
the phase shifts tend to zero, but becomes ambiguous if some of them tend 
to TT, shows that this form of the dispersive framework is deficient. As an 
illustration, we return to the simplified problem and discuss the change in the 
function Mq{s) which results if the high energy behaviour of the phase shift 
So{s) is modified. Suppose that Sq{s) is taken to rapidly drop to zero at some 
large value sq, such that the corresponding Omnes function differs from the 
original one by the factor 1 — s/sq. If the dispersion relations are written in 
the form (^, the solution merely picks up this factor and thus only receives 
a small correction at low energies. With eq.(D, however, the modification 
generates a qualitative change, as it selects a unique member from a one 
parameter family of solutions - the result cannot be trusted, because it is 
sensitive to physically irrelevant modifications of the input. 

5. Subtractions. The high energy behaviour of the amplitude is the same 
as in the case of elastic scattering and is dominated by Pomeron exchange 
(since the amplitude irrj —>■ tth violates isospin symmetry, the coupling to the 
Pomeron is suppressed by a factor of md — rriu, but this holds for the entire 
amplitude). Comparison with vrvr-scattering indicates that, in principle, the 
dispersive representation of the three functions Mo(s), Mi(s), M2(s) should 
require only two subtractions. The chiral perturbation series yields an over- 
subtracted representation: The number of subtractions grows with the order 
at which the series is evaluated. 

If only two subtractions are made, the contributions from the region above 
-ft'i^'-threshold are quite important and need to be analyzed in detail (com- 



pare e.g. ref.|]l5[, where the scalar form factors are studied by using unitarity 
also for the quasi-elastic transition utt^KK). The uncertainties associated 
with the high energy part of the dispersion integrals do not limit the accu- 
racy of the measurement of Q, however. One may use more than the minimal 
number of subtractions and determine the extra subtraction constants with 
the observed Dalitz plot distribution. With sufficiently many subtractions, 
the contributions from the high energy region become negligibly small. In 
principle, a measurement of Q requires theoretical input only for the overall 
normalization, while the remaining subtraction constants, which specify the 
dependence of the amplitude on s, t, m, may be determined experimentally. 

We do not invoke the observed Dalitz plot distribution, but use chiral 
perturbation theory to determine the subtraction constants. The number of 
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subtractions is chosen such that M (s, t, u) grows at most hncarly in all di- 
rections s,t,u ^ oo: Mo{s) = 0{s), Mi{s) = 0(l),M2(s) = 0{s). In view 
of the Omnes factors, this implies that the subtraction polynomials are 
of the form Po{s) = + PqS + •jos'^ , Pi{s) = ai+ PiS, P2{s) = a2 + P2S and 
thus involve altogether seven subtraction constants. As it is the case with 
the one loop amplitude, the decomposition into the three isospin compo- 
nents Mq{s) , Mi{s) , M2{s) is not unique. Since the transformation Mi(s) —>■ 
Mi{s) + ki, M2{s) — > M2{s) + k2-\-k^s preserves the asymptotic behaviour 
and can be absorbed in Mo(s), only four of the seven subtraction constants 
are of physical significance. Without loss of generality we exploit this free- 
dom in the isospin decomposition and set a\ = a2 = 132 = 0. The dispersion 
relations then take the form 

Mo{s) 
M2(s) 

The four subtraction constants ao? /^o, 70) A play the same role as the terms 
a, b, c, d occurring in the one loop result. 

The constants ckq, Po already appear in the current algebra approximation, 
which corresponds to ccq +/9o « = (3s — 4M^ ) / (M^ — ) . At leading order, the 
amplitude thus vanishes at s = |M^, a feature which is beautifully confirmed 
by the observed Dalitz plot distribution. The phenomenon is a consequence 
of SU(2)jjxSU(2)i-symmetry and does not rely on the expansion in powers 
of nis, which is generating the main uncertainties in the chiral representation: 
In the limit m„ = m(i = 0, where the pions are massless, the decay amplitude 
contains two Adler zeros, one at p^+=0 {s = u = 0), the other at ^,^-=0 {s = 
t = 0). The perturbations generated by m„,mrf produce a small imaginary 
part, but the line ReM{s,t,u) = still passes through the vicinity of the 
above two points. The figure shows that along the line s = u, the one loop 
corrections shift the Adler zero from 1.33 to sa = lAl M^. The position 
of the zero is related to the value of the constant ag, while (3o determines 
the slope of the curve there. The figure also shows that the corrections of 
order barely modify the latter: The slope of the one loop amplitude 
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ao + (5qS + 7oS^ H / — 

r ^\ o s ds' sin(5i(s')Mi(s' 
I 7rJ4M2 s' \Ui{s ) \ {s- s - le) 



TT J A 



's' — s — ie) 



2 ds' sin_5o(s') Mo(s') 
s ds' sin5i{s')Mi{s') 

ns -r 

2 ds' sm5'){s')M2{s') 



(10) 



4M|S'2 |02(S')| is' 



differs from the current algebra result, /?o = 3/(M^ — M^), by less than a 
percent. We therefore assume chiral perturbation theory to be reliable in the 
vicinity of the Adler zero and determine the subtraction constants and 
by matching the dispersive representation of the amplitude with the one loop 
representation there. More specifically, we require that the first two terms in 
the Taylor series 

M(s, 3so - 2s, s) = Ma+{s- sa)Sa + 0{{s - sa^} 

agree with those of the one loop amplitude, where 

Ma = -0.0041 z , 5a = (0.1957 - 0.0519 i)M-2 . 

Concerning the value of ao, this procedure is safe, because ao is protected 
by the low energy theorem mentioned above. The slope, on the other hand, is 
not controlled by SU(2)jjxSU(2)£,. The one loop result does account for the 
corrections of order rris, but disregards higher order contributions. According 
to a general rule of thumb, first order SU(3) breaking effects are typically 
of order 25% while second order contributions are of the order of the square 
of this. It so happens that in the case of the slope, the various first order 
corrections nearly cancel one another. The second order effects are discussed 



in some detail in ref.[16|, on the basis of the large expansion. In particular, 
it is shown there that the one loop result fully accounts for rji]' mixing [[17|] , 
except for a factor of cos6'^r,/ ~0.93. This represents a second order correction 
of typical size, consistent with the rule of thumb. 

In the convention adopted above, the term f3i determines the slope of 



Mi(s). More generally. Pi is related to the convention independent quantity 
M((0) + iM^'(0) = A + ds' 



1 ^ j \ sin5i(s')Mi(s') ^ sin52(s')M2(s') 



(11) 

In one loop approximation, the Ihs coincides with 6i + C2. The low energy 
theorem (|^) thus relates Pi to the coupling constant L3, whose value is known 
rather accurately from K^^ decay, = (—3.5 ± 1.1) ■ 10~^ 0. 

Both for the dispersive representation and for the one loop approxima- 
tion, the quantity M((0) + |M^'(0) receives two contributions: one from an 
integral over the discontinuities of the amplitude, which accounts for the 
effects generated by the low lying states, and one from a subtraction term 
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{(3i ^ -(4L3 - l/QAiT^)/F^{M^ - M^)}, which incorporates all other sin- 
gularities, including those at infiniy. The main difference between the two 
representations is that the integral in eq.(|TT]) exclusively accounts for the vrvr 
discontinuities, while the one in eq.(|D includes the singularities generated 
by KK, ITT] and rjr] intermediate states, albeit only to one loop. One readily 
checks that the integrals over the elastic region 4M2 < s < (M^ + only 
differ by terms of order p^, which are beyond the accuracy of the one loop 
prediction. This is an immediate consequence of the fact that chiral pertur- 
bation theory is consistent with unitarity. In the elastic region, the dispersive 
calculation yields a more accurate representation of the discontinuities than 
the one loop approximation. Accordingly, we use this approximation only to 
estimate the remainder. The one loop result shows that the contributions 
from TTi] intermediate states are proportional to and therefore tiny - in- 
elastic channels generate significant effects only for AMj^ <s<oo. The entire 
contribution from this interval to the integral in eq. amounts to a shift in 
Pi of +1.7GeV~'^, too small to stick out from the uncertainties in the term 
from L3. We conclude that jSi is dominated by this term, 

and take the above shift as an estimate for the uncertainties associated with 
the inelastic channels. 

The constant 70 is related to the convention independent combination 
M^'(O) + fiVfs (0) of Taylor coefficients and may be evaluated in the same 
manner. In particular, one may check that the elastic contributions to the 
relevant integrals coincide within the accuracy of the one loop representation 
also in this case. The only difference is that the low energy theorem (|^) 
does not contain a term analogous to L3, so that the arguments outlined 
in the preceding paragraph now imply 70 ~ 0. The uncertainties to be 
attached to this result may again be estimated with the contribution from the 
interval 4M|^ < s < 00 to the integral over the imaginary part of the one loop 
amplitude, which amounts to a shift in the value of 70 by +8.6 GeV~^. The 
net effect of the uncertainties in the subtraction constants /5i,7o is small. At 
the center of the Dalitz plot, the above shifts due to inelastic discontinuities 
increase the amplitude by about 5 %. The order of magnitude is comparable 
with the shift in /3o due to rji]' mixing, which acts in the opposite direction. 
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6. Iterative construction of the solution. Since the system is linear, the 
dependence of the solution on the subtraction constants oq, /3o, 7o and (5i 
is of the form Mi(s) = ao/i(s) + /?o5'i('S) + 7o^i(s) + jdi ki{s). The coef- 
ficients /o(s), /2(s) represent the solution for = 1, /^o = 7o = A = 
and similarly for the other terms. The solutions may be obtained itera- 
tively, starting e.g. with Mi{s) = 0. Once the iteration has been performed 
for the four fundamental solutions defined above, the values of ao, Pq for 
which the Adler condition is met are readily worked out: Denote the to- 
tal amplitudes which corresponds to the four sets of isospin coefficients by 
/(s, t, u),g{s, t, u), h{s, t, u), k{s, t, u), restrict these to the line s = u and cal- 
culate the values as well as the first derivatives at the point sa = 1.41 M^. 
The requirement that the sum M = ao f + g + h + (3ik agrees with the 
one loop amplitude there, both in value and slope, amounts to two linear 
equations for and Pq. 

The numerical results obtained with the above dispersive machinery are 
illustrated in fig.l: The solid line shows the real part of the amplitude along 
the line s = u, which connects one of the Adler zeros with the center of the 
Dalitz plot. As can be seen from this figure, the corrections to the one loop 
result are rather modest. The calculation confirms the crude estimates given 
in ref.0]. The error bars to be attached to the dispersive evaluation are 
remarkably small. They are dominated by the uncertainties in the low en- 
ergy theorems used to determine the subtraction constants. A more detailed 



account of our numerical results is in preparation . 

We are indebted to V. V. Anisovich, Hans Bebie, Jiirg Gasser and Peter 
Minkowski for useful comments and generous help with computer ethono- 
graphy and to Joachim Kambor, Christian Wiesendanger and Daniel Wyler 
for keeping us informed about closely related work fT^. 



References 

[1] D. G. Sutherland, Phys. Lett. 23 (1966) 384; 

J. S. Bell and D. G Sutherland, Nucl. Phys. B4 (1968) 315; 

For a recent analysis of the electromagnetic contributions, see 

R. Baur, J. Kambor and D. Wyler, Electromagnetic corrections to the 

decays rj Sn, preprint |hep-ph/9510396 . 
[2] J. Gasser and H. Leutwyler, Nucl. Phys. B250 (1985) 539. 



13 



[3] S.M. Roy, Phys. Lett. 36B (1971) 353; 

J. L. Basdevant, C. D. Froggatt and J. L. Petersen, Nud. Phys. B72 
(1974) 413. 

J. Stern, H. Sazdjian and N. H. Fuchs, Phys. Rev. D47 (1993) 3814. 
M. Knecht, B. Moussalam, J. Stern and N. H. Fuchs, The low energy 
TTTT amplitude to one and two loops, preprint |hep-ph /95073 1 9 . 



J. Bijnens, G. Colangelo, G. Ecker, J. Gasser and M. E. Sainio, Elastic 

TTTT scattering to two loops, preprint |hep-ph/95113971 . 

J. Bijnens, G. Colangelo and J. Gasser, Nucl. Phys. B427 (1994) 427. 

V. N. Gribov, Nucl. Phys. 5 (1958) 653; 

N. Khuri and S. Treiman, Phys. Rev. 119 (1960) 1115; 

V. V. Anisovich and A. A. Anselm, Sov. Phys. Usp. 9 (1966) 287 [Usp. 

Fyz. Nauk 88 (1966) 287]. 

A. Neveu and J. Sclierk, Ann. Phys. 57 (1970) 39. 

C. Roiesnel and T.N. Truong, Nucl. Phys. B187 (1981) 293; 

T.N. Truong, Nucl. Phys. B (Proc. Suppl.) 24A (1991) 93. 

A. V. Anisovich, Physics of Atomic Nuclei 58 (1995) 1383 [Yad. Fiz. 58 

(1995) 1467]. 

J. Kambor, C. Wiesendanger and D. Wyler, Final state interactions and 
Khuri- Treiman equations in rj ^ Sn decays, preprint |hep-ph/9509374 . 



J. F. Donoghue, On the marriage of xPT and dispersion relations, 
preprint |hep-ph/9506205 . 



B. Ananthanarayan and P. Biittiker, "Determination of the chiral pa- 
rameters ii and £2 from Roy equation analysis of mi phase shifts", in 
preparation. 

J. F. Donoghue, J. Gasser and H. Leutwyler, Nucl. Phys. B343 (1990) 
341. 

H. Leutwyler, Implications of rjr]' mixing for the decay rj Sn, [hep- 
ph/960123q . 



R. Akhoury and M. Leurer, Z. Phys. C43 (1989) 145, Phys. Lett. B220 
(1989) 258; 

A. Pich, rj-decays and chiral Lagrangians, in Proc. Workshop on rare 
decays of light mesons, Gif sur Yvette, 1990, ed. B. Mayer (Editions 
Frontieres, Paris, 1990). 
[18] A. V. Anisovich and H. Leutwyler, Measuring the quark mass ratio 
{m\ — ml) /{ml — m^) by means of rj—^Sir, in preparation. 



14 



